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,J^ ■ Abstract 

-)— > . 

2 ' In this paper, we study forward-backward doubly stochastic differential 

equations driven by Brownian motions and Poisson process (FBDSDEP in 
short). Both the probabilistic interpretation for the solutions to a class 
of quasilinear stochastic partial differential-integral equations (SPDIEs in 
^ ' short) and stochastic Hamiltonian systems arising in stochastic optimal con- 

^^ I trol problems with random jumps are treated with FBDSDEP. Under some 

O/ ' monotonicity assumptions, the existence and uniqueness results for measur- 

able solutions of FBDSDEP are established via a method of continuation. 
Furthermore, the continuity and differentiability of the solutions of FBDS- 
1^ I DEP depending on parameters is discussed. Finally, the probabilistic inter- 

pretation for the solutions to a class of quasilinear SPDIEs is given. 
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. F^. , 1 Introduction 



Nonlinear backward stochastic differential equations with Brownian motion as 
noise sources (BSDEs in short) were first introduced by Pardoux and Peng [lOj . 
By virtue of BSDEs, Peng [14j has given a probabilistic interpretation (nonlinear 
Feynman-Kac formula) for the solutions of semilinear parabolic partial differential 
equations (PDEs in short), for more detailed, the reader is referred to Darling and 
Pardoux [5], Pardoux and Zhang [13] and so on. Fully coupled forward-backward 
stochastic differential equations (FBSDEs in short) can provide a probabilistic 
interpretation for the solutions to a class of quasilinear parabolic and elliptic 
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PDEs (cf Pardoux and Tang |12j l and have been investigated deeply. FBSDEs 
were studied first by Antonelli [1] and Ma et al. [9] , Cvitanic and Ma [1] and used 
to hedge options involved in a large investor in financial market. Consequently, 
FBSDEs were developed in Hu and Peng [7j, Peng and Wu [15j and Peng and 
Shi [IB] and so on. 

A class of backward doubly stochastic differential equations (BDSDEs in 
short) was introduced by Pardoux and Peng [TT| in 1994, in order to provide 
a probabilistic interpretation for the solutions to a class of quasilinear stochastic 
partial differential equations (SPDEs in short). Due to their important signifi- 
cance to SPDEs , the researches for BDSDEs have been in the ascendant (cf. Bally 
and Matoussi [2], Zhang and Zhao |;28j, Ren et al. [18], Hu and Ren j6j and their 
references). In 2003, Peng and Shi [ITj have introduced a type of time-symmetric 
forward-backward stochastic differential equations, which is a generalization of 
stochastic Hamilton system. Recently Zhu et al. [29] have extended the results 
in [IT] to different dimensional forward-backward doubly stochastic differential 
equations (FBDSDEs in short) and weakened the monotone assumptions. How- 
ever, the theory of FBDSDEs has not been investigated enough so far. 

BSDEs driven by Brownian motions and Poisson process (BSDEP in short) 
was first discussed by Tang and Li [2T]. After then Situ [19j obtained an existence 
and uniqueness result with non-Lipschitz coefficients for BSDEP, so as to get the 
probabilistic interpretation for solutions of partial differential-integral equations 
(PDIEs in short). Using this kind of BSDEP Barles et al. |3j and Yin and Mao 
|25j proved that there exists a unique viscosity solution for a system of parabolic 
integral-partial differential equations. Fully coupled forward-backward stochastic 
differential equations with Poisson process (FBSDEP in short) was discussed by 
Wu [22] and Yin and Situ [26]. Then FBSDEP were used to study the linear 
quadratic optimal control problems with random jump by Wu and Wang [23] 
and the nonzero-sum differential games with random jumps by Wu and Yu [24]. 
Recently BDSDEs driven by Brownian motions and Poisson process (BDSDEP 
in short) was discussed by Sun and Lu [20] . 

Because of their important significance to SPDEs, it is necessary to give in- 
tensive investigation to the theory of FBDSDEs. In this paper we study FBDS- 
DEs driven by Brownian motions and Poisson process (FBDSDEP in short), 
which generalize the so-called time-symmetric forward-backward stochastic dif- 
ferential equations introduced by Peng and Shi [T7]. FBDSDEP can provide more 
extensive frameworks for the probabilistic interpretation (nonlinear stochastic 
Feynman-Kac formula) for the solutions to a class of quasilinear stochastic par- 
tial differential-integral equations (SPDIEs in short) and stochastic Hamiltonian 
systems arising in stochastic optimal control problems with random jumps. Un- 
der some monotonicity assumptions, we establish the existence and uniqueness 
results for measurable solutions of FBDSDEP by means of a method of continu- 
ation systemically introduced by Yong [27J. Then we discuss the continuity and 
differentiability of the solutions to FBDSDEP depending on parameters. Fur- 
thermore, by virtue of FBDSDEP, we give the probabilistic interpretation for the 
solutions to a class of quasilinear SPDIEs. Finally, we discuss a doubly stochastic 
Hamiltonian system. 



The paper is organized as follows. In Section [21 some preliminary results are 
given. Section [3] is devoted to proving the existence and uniqueness result for 
FBDSDEP. In Section U the continuity and differentiability of the solutions to 
FBDSDEP depending on parameters is discussed. In Section [Sj the probabilistic 
interpretation for the solutions to a class of quasilinear SPDIEs is given by virtue 
of this class of FBDSDEP. Finally, in Section [U the above results are applied to 
a doubly stochastic Hamiltonian system. 

2 Preliminary 

Let {0,,T,P) be a complete probability space, and [0,T] be a fixed arbitrarily 
large time duration throughout this paper. We suppose {J-'t}t>Q is generated by 
the following three mutually independent processes: 

(i) Let {Wi;0 <t < T} and {Bt;0 <t< T} be two standard Brownian mo- 
tions defined on (0, J-", P), with values respectively in M and in M . 

(ii) Let A^ be a Poisson random measure, on M_|_ x Z, where Z C M*" is 
a nonempty open set equipped with its Borel field B{Z), with compensator 
N{dzdt) = X{dz)dt, such that N{A x [0,t]) = (N - N){A x [0,t])t>o is a martin- 
gale for all A € B{Z) satisfying X{A) < oo. A is assumed to be a cr-finite measure 
on {Z,B{Z)) and is called the characteristic measure. 

Let Af denote the class of P-null elements of T. For each t G [0, T], we define 
J^t = J^t^ V T^rp V J^/^, where for any process {rjt}, J^^^ = a {r]r — Vs] s < r < t} V 
Af,J^ = J'of Note that the collection {Tt,t E [0,T]} is neither increasing nor 
decreasing, and it does not constitute a classical filtration. 

We introduce the following notations: 

M2(0,T;M") = {t;i,0<t<T, is an M'^-valued, Ji-measur able process 

fT 
such that E \vt\ dt < oo}, 
Jo 

F;v(0,T;M") = {kt,0 < t <T, is an M"- valued, Jt-measurable process 

such that E / |A;j(2)| X{dz)dt < oo}, 

Jo J z 

Ly^,\{MP') = {k{z), k{z) is an M"- valued, ^(Z)-measurable function 

such that ||A;|1 = ( / \k{z)fX{dz))'^/'^ < oo}, 
Jz 

L (rj, J^T, P;]R"') = {^, .^ is an M"-valued, J-^'-measurable randoin 

variable such that E\S^\^ < oo}. 



For a given u € M'^{0, T; R'^) and v G M'^{0, T; M'), one can define the (standard) 
forward Ito's integral jQUgdWg and the backward Ito's integral / VgdBs- They 
are both in M^(0,T;M) (see [H]). We use the usual inner product (•,•) and 
Euclidean norm | ■ | in M", M"^, M"^^' and R"^"^. Ah the equalities and inequalities 
mentioned in this paper are in the sense of dt x dP almost surely on [0,T] x Q.. 



Consider the following BDSDE with Brownian motions and Poisson Process 
(BDSDEP in short): 

' dPt = F{t,Pt,Qt,Kt)dt + Git,Pt,Qt,Kt)dBt 

+QtdWt+ f Kt_iz)N{dzdt), (1) 

Jz 
Pt = e, 



where 



F : [7 X [0,T] X M™ X R""^^ x L^(.)(M'") 

-A(.)( 



G : J7 X [0, T] X M™ X R""^^ x h] ^™™^ ^ ™™^' 



2^n 7r„ P;M™); 



Definition 2.1 >1 triple of J^t-'^n.easurahle stochastic processes {P,Q,K) is called 
a solution to BDSDEP ^^ if {P, Q,K) G M^{[0,T];R"'+"''"^) x F^{0,T;M."') 
and satisfies BDSDEP |ip. 

We assume that 

(HI) {^ CeL^{n,TT,P 

(ii) F is J^i-progressively measurable and satisfies F(a;, t, 0, 0, 0) € M^(0, 

r;M"^); 

(iii) G is 7"t-progressively measurable and satisfies G(a;, t,0, 0, 0) G M^(0, 

(iv) F and G satisfy Lipschitz conditions to P, Q, K, that is, there exist 
c > and < 7 < 1 such that 

\F{t,Pi,Qi,K,) - F{t,P2,Q2,K2)f 

< c{\Pi - P2P + \Qi - Q2I' + lli^i - K2f), 
\Git,Pi,Qi,Ki)-G{t,P2,Q2,K2)\^ 

< c\Pi - P2I' + 7(IQi - Q2I' + lli^i - K2f)- 

In order to attain our results, we give the following Proposition 12 . 21 and Propo- 
sition [2]3] appeared in pO] and Proposition 12.41 appeared in [19]. 

Proposition 2.2 Under the assumption (HI), BDSDEP ([I]) has a unique solu- 
tion {P,Q,K) G M2([0,r];M™+™>^'=') x F^{0,T;R"'). 

Proposition 2.3 Let a G A'P{[0,T];R"'), /3 G M'^{0,T;R'^),j G M'^{0,T; 
M™^'). 6 G M2(0,r;IR""^'^) be such that: 

at = ao + / l^sds + / 7sdSs + / SgdWs 
Jo Jo Jo 



+ Ks^{z)N{dzds), 0<t<T. 

Jo Jz 



Then 

ft rt 

|2 |„ |2 



\at\ = |ao| +2/ {as,Ps)ds + 2 I {as,jsdBs) 
Jo Jo 

+2 I {as,5sdWs) + 2 I [ {a,,Ks_{z))N{dzds) 
Jo Jo Jz 

- / \-fs\^ds+ / \6s\^ds+ / \\Ksfds, 
Jo Jo Jo 

E\at\'^ = E\ao\'^ + 2E I {as,l3s)ds-E I \-fs\'^ds 

Jo Jo 

+E [ \5s\'^ds + E [ \\Ksfds. 
Jo Jo 

Proposition 2.4 Let a G M2([0,T];IR"),/3 G M^{0,T;W'),j G M'^{0,T;W'^] 
S e F^(0,r;M") 6e suc/i i/iai; 

at = ao+ I Ms + I TsdWs + 6s_{z)N{dzds), 0<t<T. 

Jo Jo Jo Jz 

Then for all u{s,x) G C^'^{[0,T] x M";M'"), one has that 

u{t,at) - u{0,ao) = / Cu{s,as)ds+ / Vu{s,as)jsdWs 

Jo Jo 

+ / / {u{s-,as_ + 6{s-,z)) - u{s-,as_))N{dzds), 
Jo Jz 



where 



Lui 
Cu 

LUr. 



with 

Luk{t,x 
duk 



(t,x) + j;ft(t)|^(t,x) + i Y. (^^*hit)i:^it, 



X 



dt ^ dxi 2 .■^-^^ dxidx^ 

/* ^ ^ 

+ / (iifc(t,x + (5(t,x,z)) -Uk{t,x) -'YSi{t,x,z)——{t,x))X{dz), 
Jz ^^i Xi 

k = I,- • • ,m. 



3 The existence and uniqueness theorem of FBDS- 
DEP 

Consider the following forward-backward doubly stochastic differential equations 
with Brownian motions and Poisson process (FBDSDEP in short): 

' dXt = fit, Xt, Pt, Yt, Qt, Kt)dt + g{t, Xt, Pt, Yt, Qt, Kt)dWt - YtdBt 
h{t.,Xt_,Pt^,Yt_,Qt_,Kt_iz),z)N{dzdt), 



dPt = F{t, Xt,Pt, Yt, Qt, Kt)dt + G{t, Xt,Pt, Yt, Qt, Kt)dBt + QtdWt (2) 



+ / Kt_{z)N{dzdt), 
Xo = ^{Po), Pt = HXt) 



where 
F 

f 
G 

g 

h 



n X [0, T] X M" X M™ X 
n X [0, T] X M" X M™ X 
n X [0, T] X M" X M™ X 
n X [0, T] X M" X M" X 
n X [0, T] X M" X M™ X 



tnxl 
hnxl 
tnxl 
pnxl 









*) X Z 



r, 

^mxl 
nnxd 



r, $ : SI X M" ^ u. . 
Given an ?tt, x n full-rank matrix H. Let us introduce some notations 

U = {X, P, Y, Q, K), A{t, U) = {H^F, Hf, H^G, Hg, Hh){t, U), 
{A,U) = {X,H^F) + {P,Hf) + {Y,H^G) + {Q,Hg) + {{K,Hh)), 

where 

H^G = {H^Gi---H^Gi), Hg = {Hgi---Hgd), 

{{K,Hh)) = [ {Kt{z),Hh{t,Ut,z))X{dz). 



Definition 3.1 ^ quintuple of Ft-measurable stochastic processes {X, P, Y, Q, K) 
is called a solution to FBDSDEP (0J, if {X,P,Y, Q,K) G M'^{Q,T; 
^n+m+nxi+mxd^^ ^ F^(0,r;M™) and satisfies FBDSDEP (0). 

The following monotonicity conditions are our main assumptions: 

(H2) VC/ = {X,P,Y,Q,K),U = {X,P,Y,Q,K) G M"+™+'^>^'+™><'^ x L\^^.^{W^), 
Vte[0,T], 

{A{t,U)-A{t,U),U-U) 

< -^,,(\H{X-X)\' + \HiY-Y)\') 

-f^2i\H^iP -P)\' + \H^iQ - Q)P + 11^^ (^ - K)f)- 



(H3) (^(P) -*(P),iJ'^(P-P)) < -/32|i?'^(P-P)P, VP,PgM'", 

($(x) -^{x),H{x -X)) >i3i\H{x-x)\'^, yx,x gM". 

Here, /xi, fi2, Pi and /32 are given nonnegative constants with /ii + /i2 > 0, 
/3i + /?2 > 0, /xi + /32 > 0, /U2 + /?! > 0. Moreover we have /xi > 0, /3i > (resp. 
/X2 > 0, /32 > ) when m > n (resp. n> m). 
We also assume that 

(H4) For each [/ G M«+'"+"x'+™x«' x L^^ ^(M™), A(-, ^7) is an Ji-measurable vec- 
tor process defined on [0,r] with A{-,0) G M^{0,T; M.'^+m+nxi+mxd^ ^ 
F^(0, T; M™-), and for each X G M", ^(X) is an Jr-measurable vector pro- 
cess with $(0) G L2(0,7't,-P;M™), and for each P G M"", *(P) is an 
Jo-measurable vector process with ^(0) G L'^{i^,To,P;Mr'). 

(H5) /, F,g, G, h, ^ and ^ satisfy the Lipschitz conditions: there exist constants 
c> and < 7 < 1 such that \/U = {X, P, Y, Q, K), U = {t, X, P, Y, Q, K) 

^ ^n+m+nxl+mxd ^ ^^^.^ (R'"), Vt G [0,r], 

\f{t,X,P,Y,Q,K)-f{t,X,P,Y,Q,K)\^ 

< c{\x - x\^ + \p- p\^ + \y - y?' + \Q- Q? + \\K - Kf), 

\F{t, X, P, Y, Q, K) - F{t, X,P,Y,Q,K)\^ 

< c{\X - Xp + \P- Pp + 1^ - ?P + \Q- OP + 11^ - Kf), 
\g{t,X,P,Y,Q,K) - g{t,X,P,Y,Q,K)f 

< c{\x - x|2 + \P- -Pp + \Q- Qp + \\K - Kf) + 7|y - yp, 

\G{t,X,P,Y,Q,K)-G{t,X,P,Y,Q,K)\^ 

< c{\X - X|2 + \P- P|2 + \Y- y|2) + 7(|g - Q\^ + \\K - Kf), 
\h{t,X,P,Y,Q,K)-h{t,X,P,Y,Q,K)\^ 

< c{\X - Xp + \P- Pp + \Y - yp + \Q- Qp + \\K - Kf), 
|^(P) - ^(P)| < c\P - P\, \^{X) - ^{X)\ < c\X - X\. 

Then we claim the main result of this section. 

Theorem 3.2 Under the assumptions (H2)—(H5), ^ has a unique adapted so- 
lution {Xt,Pt, Yt,Qt,Kt) G M^{0,T; ]Rn+™+"x/+mxd) ^ F^{0,T;R"'). 

The proof of this theorem is divided into two parts, i.e., existence and uniqueness. 
We first give the proof of uniqueness. 

Proof. (Uniqueness) Let Us = {Xs,Ps,Ys,Qs,Ks} and U^ = {X'„P;„Y^,Q'^,K'^) 
be two solutions of ©. We set U = {X,P,Y,Q,k) = {X - X',P - P',Y - 



Y',Q- Q',K - K'). Applying Ito's formula to {HXs,Ps) on [0,T], we have 
E{HXt, HXt) - HX!r)) - E{H{^{Po) - ^(P,^)), Pq) 

= E [ {A{s,Us)-A{s,U',),Us)ds 
Jo 

< -niE [ {\HX,\'^ + \HY,\'^)ds 
Jo 

-fi2E / (|if^P,|2 + \H^Qs\'^ + \\H'^Ksf)ds. 
Jo 



Then 



T 




l^iE I {\HXs\^ + \HYs\^)ds 



+fi2E / {\H^Ps\'^ + \H'^Qs\'^ + \\H'^Ksf)ds < 0. 
Jo 

If m > n and fii > 0, then we have \HXt\'^ = 0, \HYt\'^ = 0. Thus Xt = 
X^,Yt = y/. In particular, ^{Xt) = ^{X!p). Consequently, from Proposition 12.21 
it follows that Pt = P/, Qt = Q't and Kt = K[. 

li m < n and //2 > 0, then we have Pt = Pt,Qt = Q't-, and Kt = K[. In 
particular, ^(Pq) = *(Po)- Thus from Proposition [221 it follows that Xt = X[, 
and Yt = F/. 

Similarly to above arguments, the desired result can be obtained easily in the 
case m = n. D 

Remark 3.3 In the proof of the uniqueness and existence, (H2) and (H3) can 
he replaced by 

(H2)' VC/ = {X,P,Y,Q,K), tJ = {X,P,Y,Q,K) G wi+m+nxi+mxd ^ /.^^.^(M™), 
VtG[0,T], 

{A{t,U)-A{t,U),U-U) 
> fii{\H{X-X)\^ + \HiY-Y)\^) 

+li2{\H^{P -P)? + \H^(Q -Q)? + \\H^{K - K)f), 

(H3)' (*(P) -*(P),P''^(P-P)) >/32|F'^(P-P)P, yP,P(£R'^, 

{<^{x) -^x),H{x -X)) < -i3i\H{x -x)]"^, yx,x gM". 

Here, ^i, ^2, /3i and /32 ore given nonnegative constants with /Ui + ^2 > 0, 
/3i + /32 > 0, /^i + /32 > 0, /U2 + /3i > 0. Moreover we have fii > 0, f3i > (resp. 
/i2 > 0, /32 > J when m > n (resp. n > m). 

The proof of the existence is a combination of the above technique and the 
method of continuation systemically introduced by Yong [27]. We divide the 
proof of existence into three cases: m > n, m < n and m = n. 

Case 1 m > n 



If m > n, then ^ui > and /3i > 0. We consider the following family of 
FBDSDEP parametrized by a G [0, 1]: 

' dXt = [a/(t, Ut) + /o(t)] dt + [ag{t, Ut) +go{t)\ dWt - YtdBt 
+ J^[ah{t^,Ut_,z) + ho{t^,z)] N{dzdt) , 
dPt = [aF{t,Ut)-{l-a)fiiHXt + Fo{t)]dt + QtdWt ^ (3) 

+ [aG{t, Ut) - (1 - a)fiiHYt + Go{t)] dBt + f^ Kt_ {z)N{dzdt), 
Xo = a^(Po) + V-, Pt = a^Xr) + (1 - a)HXT + (f, 

where Ut = {Xt, Pt,Yt,Qt), [Fo, fo,Go,go,ho) G Af2(0, r;M-+"+™x'+"x^) x 
F|r(0,T;M"), V' G L'^{n,To, P;W) and 93 G ^^(Jl, J-^, P;M™) are given arbi- 
trarily. 

When a = 1 the existence of the solution of ([3]) implies clearly that of ([2]). 
Due to Proposition 12. 2| when a = 0, the equation ([3]) is uniquely solvable. The 
following apriori lemma is a key step in the proof of the method of continuation. 
It shows that for a fixed a = oq € [0,1), if ([3]) is uniquely solvable, then it is 
also uniquely solvable for any a G [ao,ao + ^0]) for some positive constant 60 
independent of oq. 

Lemma 3.4 We assume m > n. Under the assumptions (H2)-(H5), there exists 
a positive constant 5q such that if, apriori, for each ip G L^(ri, J-q, P;M"), if G 
L^in,TT,P;^n, and (Fo,/o,Go,5o,/io) S M2(0,T;R-+"+™x'+"X'^) x F|,(0,T; 
M"^), (^ is uniquely solvable for some ao G [0, 1), then for each a G [ao, ao + 5o], 
and V G L2(J], Jo,P;M"), ^ G L^n,TT,P;Rn, {Fo, fo,Go,go,ho) G M^{0,T; 
^m+n+mxi+nxd^^ X F^{0,T;W^), ^ is also uniquely solvable in M^{0,T; 

^n+m+nxl+mxd^ ^ F^(0, T; M""). 

Proof. Since, for each ^ G L'^{n,Fo,P;W'), 93 G ^^(f], J-^., P;M'»), (Fo,/o, 
Go,go,ho) G M2(0,r;M™+'^+™x'+"xd) ^ F^{0,T;R"') there e_xists a unique so- 
lution of dSD for a = ao. Thus, for each U = {X,P,Y,Q,K) G M'^{0,T; 
^n+m+nxi+mxd^ xF^(0, T; M*"), there exists a unique U = {X,P,Y,Q,K) G 
M2(0,T;M"+"+"x'+™xc() X F2,(0,T;R'") satisfying the following FBDSDEP: 

dXt = [aof{t, Ut) + 6 fit, Ut) + /o(t)] dt + [aog{t, Ut) + Sg{t, Ut) + go{t)] dWt 

-YtdBt + / [aoh{t^,Ut_,z) + 6h{t^,Ut_,z) + ho{t^,z)]N{dzdt), 
Jz 

dPt = [aoF{t,Ut)-{l-ao)fiiHXt + d{F{t,Ut)+fiiHXt) + Fo{t)]dt 

+ [aoGit, Ut) - (1 - ao)fiiHYt + 5(G(t, Ut) + fnHYt) + Go(t)] d^t 

+QidWt+ f Kt_{z)N{dzdt), 
Jz 

Xo = ao^iPo) + 6^iPo)+i^, 

Pt = ao^{XT) + {l-ao)HXT + 6{^{XT)-HXT) + ip, 



where 6 € (0, 1) is independent of a^. We will prove that the mapping defined by 

{Ut,XT,Po)=Iao+5{Ut,XT,Po): 

M2(0,r;IR"+™+"^'+™^'^) X F^(0,r;IR™) 
^ M2(0,r;E"+™+"^'+™^'^) X F^(0,r;E™) 

is contractive for (5 > which is small enough. 
Let 

(x;, p;, y/, Qj, KO = [/; = /,„+5(f7;), 

= (Xi - Xl Pt - P;, >"* - Y^, Qt - Qt, Kt - K^), 

Ut = Ut-Ul = {Xt,Pt,Yt,Qt,Kt) 

= {Xt - Xl Pt - P;, Yt - Yi Qt - Q't, Kt - K't). 

Applying Ito's formula to {HX,P) on [0, T] yields 

E{HXt, ao$(^T) + (1 - ao)HXT) 

-E [ {ao{A{s,Us)-A{s,Us)),Us)ds 
Jo 



where 



+{l-ao)fiiE [ {\HXs\^ + \HYs\^)ds 
Jo 

E{HXt, 5HXt) - E{HXt, 6^{Xt)) + E{H{ao^{Po) + 6^{Po)), Po) 

+SE f ((P„ Hf) + {HXs,F) + {HYs, G) + (4, Hg) + ((K„ Hh)))ds 
Jo 

+SiiiE [ {{HXs,HXt) + {HYs,HYs))ds, 
Jo 



f = fit, Ut) - fit, Ui), g = git, Ut, Kt) - git, U^), 
F = Fit, Ut) - Fit, tj't), G = Git, Ut) - Git, tj't), 
h = hit,Ut,-)-hit,Ul,-), 

$(Po) = ^(Po) - ^iP^), $(Po) = ^(Po) - ^iP^), 
$(Xt) = $(Xt) - $(X^), $(Xt) = $(Xt) - ^iX^) 



10 



Let B^ = m2(0,T;M"+™+"^'+™^'^) x F^(0,r;IR"^). Noting m > n, by virtue of 
(H2)-(H5), we can easily deduce 



(1 - ao + aoPi)E\HXT\'^ + f^E I {\HXs? + \HYs\^)ds 

Jo 

< 6CE I {\\Us\\l2 + \\Us\\l2)ds 
Jo 

+5C{E |PoP + E\Xt\^ + E\Pq\^ + E\Xt\^). (4) 

with some constant C > 0. Hereafter, C will be some generic constant, which 
can be different from line to line and depends only on the Lipschitz constants c, 
7, yUi, /3i, H and T. It is obvious that 1 — ao + ao/3i > /3, /3 = min(l, /3i) > 0. 

On the other hand, for the difference of the solutions (P, Q, K) = (P — P', Q — 
Q',K — K'), we apply a standard method of estimation. Applying Proposition 
T3to |Pt|2 on [t,T], we get 

E\Pt\^ + E \Q,\^ds + E WKsfds 

= E\aQ^{XT) + (1 - ao)HXT + 5{^{Xt) - HXt)? 

-2E I {Ps,aoF{s, Us) - (1 - ao)fiiHX, + 5{F{s, Us) + iJLiHXs))ds 

+E I |aoG(s, [/,) - (1 - aQ)fiiHYs + 5{G{s, Us) + fiiHYs)\^ds. 

By virtue of (H5), we have 

1 — f^ 

E\Pt\^ + ^E {\Qs\^ + \\Ksf)ds 

< CE \Psfds + C{E\XT\'^ + 6E\XT\'^) 
+CE I {\Xs\^ + \Ys\^ + 6\\UsfB2)ds. 



It 
Then we can deduce 

f-T 



E\Po\^ + E [ {\p^\^ + \Q^\^ + \\Ksf)ds 
Jo 

l-T ^ ^ ^ 

< C{E\XT\^ + 6E\XTf) + CE {\Xs\^ + \Ys\^ + 6\\Us\\l2)ds. (5) 

Jo 

Combining the above two estimates ([4]) and ([5]), for a sufficiently large constant 
C > we can easily have 



E [ \\Us\\l2ds + E\Xt\^ + E\Po\^ 
Jo 

< 5C{E f \\Us\\l2ds + E\XT\^ + E\PQ\^), 
Jo 
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We now choose 60 = — ;. It is clear that, for each fixed 6 G [0, 5o]i the mapping 
lao+s is contract in the sense that 

E \\Us\\l2ds + E\Xt\^ + E\Po\^ 
Jo 

1 r^ -^ 

< o(^ / ||t/.|||2C^S + ^|^T|' + ^|^0|'), 

^ Jo 

Thus, this mapping has a unique fixed point U = {X,P,Y,Q,K) € M^(0,T; 
^n+m+nxi+mxd^ xF^(0,r;M™), which is the solution of ([3]) for q = ckq + S, 

Se[o,6o]. a 

Case 2 m < n 

li m < n, then ^2 > and /32 > 0. We consider following equations: 

dXt = [af{t, Ut) - (1 - a)fi2H^Pt + /o(t)] dt 

+ [a5(i, Ut) - (1 - a)ii2H^Qt + 5o(t)] 'iVFt - Ytd^i 

+ / [ah{t^,Ut_,z) + (1 - a)/i2^'^i^t_(^) + /io(i-,^)] N{dzdt), 

dPt = [aPtt, Ut) + Fo(t)] dt + QtdWt + [aG(t, Ut) + Go(t)] dS* 

- / Kt_iz)N{dzdt), 

Xo = a^{Po) + (1 - a)H^Po + V', -Pr = a$(XT) + if. 

(6) 
Due to Proposition 12.21 when a = 0, the equation ([6]) is uniquely solvable. 
When a = 1 the existence of the solution of ([6]) implies clearly that of ^. By 
the techniques similar to Lemma 13.4^ We can also prove the following lemma. 

Lemma 3.5 We assume m < n. Under assumptions (H2)-(H5), there exists 
a positive constant 5q such that if, apriori, for each ip G L^(il, J-o,P; i?"), ip G 
L^n,TT,P;Rn, and (Fo,/o, Go, 50, Me^'(0,T;i?"+"+-x'+"X'^)xF|.(0,T; 
M"^), ^ is uniquely solvable for some oq G [0, 1), then for each a G [ao, ao + 5o], 
and i^ G L2(f),7-o,P;i?"), ^ G L^{n,TT,P;R^), {Fo, fo,Go,go,ho) G M^{0,T; 
Rm+n+mxi+nxd-^ X F^ {0 , T ; R"") , ^is also uniquely solvable in M^{0,T; 

j^n+m+nxl+rnxd^ ^ F^(0,r;M™). 

Case 3 m = n 

From (H2) and (H3), we note that we only need to consider two cases as 
follows: 

(1) If /Ui > 0, /i2 > 0, /3i > 0, and /32 > 0, we can have the same result as 
Lemma 13.41 

(2) If fii > 0, ^2 > 0, (3i > 0, and /32 > 0, the same result as Lemma [37 
holds. 

Now we can give proof of the existence of Theorem 

Proof of the existence of Theorem 13.21 For Case 1, we know that, for 

each ^ G L^{n,J^o,P;Rn, V e L^{n,TT,P;Rn, {FoJo,Go,go,ho) G M^{0,T; 
j^m+n+mxi+nxd^ ^ ^^.(O, T; M™), ([3]) lias a unique solution as a = 0. It follows 
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from Lemma 13.41 that there exists a positive constant 5q = (5o(c, 7, /3i,;Ui, H,T) 
such that for any 6 G [0,6o] and ip G L'^{n,To,P;R'), 99 G L'^{il,TT,P] R""), 
{FoJo,Go,go,ho) G M2(0,r;i?™+"+'">^'+"^'^) x F^{0,T;R"'), © has a unique 
solution for a = 6. Since (5o depends only on c, 7, /3i, ^1, ff, and T, we can 
repeat this process for A^ times with 1 < A^^o < 1 + 60. In particular, for 
a = 1 with {Fq, fQ,Go,go,hQ) = 0, 99 = 0, ^ = 0, ([3]) has a unique solution in 

For Case 2, we know that, for each ip G L^(il,7o,P;M"'), 99 G L'^{Q,,Ft, 
P; M™), (Fo,/o,Go,5o,/io) G M\0,T; Rm+n+mxi+nxd^ X F2,(0,r;M™), ([6]) 
has a unique solution as a = 0. It follows from Lemma 13.51 that there ex- 
ists a positive constant 60 = 6o{c,'y, P2, IJ'2,H,T) such that for any 6 G [0,5o] 
and V G L\n,Fo,P;R^), ip G L2(Q, J-r,P; i?™), (Fq, /o, Go,5o, /io) G M2(0,r; 
^m+„+mx«+nxd) X ^2^^Q^ j^. j^m)^ p j^^^g ^^ unique solution for a = 6. Since do de- 
pends only on c, 7, /32, ^2; ^) and T, we can repeat this process for A^ times with 
1 < N60 < 1 + (5o- In particular, for a = 1 with (Fq, /o, Gc^O; ^0) = 0, 99 = 0, 
V' = 0, ® has a unique solution in M2(0,r; i?"+™+'^>^'+™><'^) xF2,(0,T;R'"). 

Similar to the above two cases, the desired result can be obtained in Case 3. 

n 

4 FBDSDEP depending on parameters 

In this section, the continuity of the solutions to FBDSDEP depending on pa- 
rameters is discussed. Let {fa, da, ha,Fa,Ga, ^a; <l*o, a G i?} be a family of data 
of FBDSDEP as follows 

dXf = Ut, X^,Pt^,Yt^,Qf,Kndt + gait, Xf , Ff , y,", Qf, K'^)dWt 
-Y,-dBt + J K{t.,Xf^ , F° , y," , Qf_ , Kl (z), z)Nidzdt), 

dPt" = Fait, Xf , FfVt", Qf,Kndt + Gait, Xf , Ff , y,", Qf, K^)dBt (7) 
+QtWt+ f KtJz)Nidzdt), 

with solutions denoted by (X°,F",y°,Q",F:'^). 
Let us give some assumptions: 

(H6) (i) The family {fa,Fa,ga, Ga, ha, ^a, ^a',Oi £ R} satisfy the equi- 
Lipschitz conditions with the same constant c as in (H5); 

(ii) {fa,Fa,ga, Ga, ha, ^a, ^a} are continuous for a, in their existing 
space norm sense respectively. 

Then we have the following continuity result. 

Theorem 4.1 Let {fa, Fa,ga,Ga, ha, ^a, ^a', a G M} 6e a family of FBDSDEP 
satisfying (H2)- (H6) with solutions denoted by (X", F°, y°, Q°, i^°) . Then, 
the family of functions {(X", F",y", (5",i^",X^,F^); a G R} is continuous 
for a in M^iO,T; ]^n+m+nxi+mxd-^ ^ F^(0,r;M™) x L2(0,Ft,F;]R") x L'^iQ, 
Fr,F;R™) . 
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Proof. For notational convenience, we only prove the continuity of FBDSDEP 
dZD at a = 0. We need to get that (X",P",y", Q",i^",X|, P(f ) converges 
to (XO,pO, yO,QO,KO,XO,pO) in M2(0,r;R"+™+"^'+'">^'^) x F^{0,T;W^) x 
L2(ri,jr^,P;R") X L2(fi, JTy, P;R'") as a ^ 0. 

We set 

Ut = Ut"-Uf = {Xt,Pt,Yt,Zt,Kt) 

= {X- - X?, p- - po, y,° - y,o, Qf - Ql K- - K?). 

Thus 

r dXt = [fait, Un - Mt, U^))dt + {g^{t, Un - go{t, U?))dWt - YtdBt 



+ 1 {Kit.,Ut"_,z) - hoit.,Ul'_z))N{dzdt), 
dPt = (F„(t, Un - Fo{t, Uf))dt + (G„(t, C/f ) - Go(t, U^))dBt 



+QtdWt + / K^ iz)N{dzdt), 
Jz 



[ Xo = ^a{P^) - ^o{PS), Pt = ^a(X^) - MX!^)- 

Applying Ito's formula to iX^p, jPtP, and (HXt,Pt) with the usual technique 
similarly to Lemma [231 we can obtain 

/■^ ^ 

E WUsWlids + EIXtI"^ + E\Po\'^ 
Jo 

< CE f {\Us,U^)-Us,U^)\^ + \gUs,U^)-gois,U^)\^)ds 
Jo 

+CE [ i\F^is,U^)-Fo{s,U^)\^ + \GUs,U^)-Gois,U^)\^)ds 
Jo 

+CE f \\ha{s,Ul-)-ho{s,Us,-)fds 
Jo 

+CE|$„(yO ) - $o(y^)P + CEl^^iY,^) - *o(>^°)P. 

where C > is some constant depending on the Lipschitz constants c, 7, fii, 
H2, /3i, /32, H, and T. Therefore, from (H6), it follows that (X",P°,y°,Q", 
K'^,X!^,P^) converges to {X^, P^,Y^,Q^,K^,X^, P^) in the space M^{0,T; 

j^n+m+nxl+mxd^ ^ ^2^(q^ rp. jgm) ^ ^2^^^ j-^^ p. pn^ ^ j2^^^ -p^^ p. pm^ aS a ^ 0. 

n 

In fact, we can also discuss the differentiability of the solution to FBDSDEP 
depending on parameters. The method is similar. These two properties are 
important and make FBDSDEP be used widely especially in practice. 

5 The probabilistic interpretation of SPDIEs 

The connection of BDSDEs and systems of second-order quasilinear SPDEs was 
observed by Pardoux and Peng |llj. This can be regarded as a stochastic version 
of the well-known Feynman-Kac formula which gives a probabilistic interpretation 
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for second-order SPDEs of parabolic types. Thereafter this subject has attracted 
many mathematicians, referred to Bally and Matoussi [2], Zhang and Zhao [28], 
Hu and Ren [6], see also Ren et al. [S]]. In [IHl]i the authors got a probabilistic 
interpretation for the solution of a semilinear SPDIE, via BDSDEs with Levy 
process. This section can be viewed as a continuation of such a theme, and 
will exploit the above theory of fully coupled FBDSDEP in order to provide a 
probabilistic formula for the solution of a quasilinear SPDIE. 

For each x G M", let {Xt,Pt,Qt,Kt;0 <t<T} denote the solution of the 
FBDSDEP: 



dXt 
dPt 

. Xt 
where 



/(t, Xu Pu Qt, Kt)dt + g{t, Xt,Pt)dWt + f h{t^,Xt_ , z)N{dzdt), 
F{t, Xt, Pt, Qt, Kt)dt + G{t, Xt,Pt, Qt, KtfdBt 



+QtdWt + / Kt_ 
X, Pt = ^{Xt), 



{z)N{dzdt), 



(8) 



F 

f 
G 

9 
h 



[o,r]x 
[o,r]x 

[0,T]x 

[o,r]x 

[0,T]x 



X 
X 



r X z 






pmxl 
hnxd 
pn 



satisfy (H2)-(H5), and 

(H7) /, g, h, F and G are of class C^, and $ is of class G^. 

We now relate FBDSDEP ([8|) to the following system of quasilinear second- 
order parabolic SPDIE: 

Cu(t,x)dt 

= F(t, X, u{t, x), Vu{t, x)g{t, X, u{t, x)),u{t, x + h{t, x, •)) — u{t, x))dt 
+G{t, x, u{t, x), Vu{t, x)g{t, X, u{t, x)),u{t, x + h{t, x, •)) — u{t, x))dBt, 
^ (t,x) G [0,r] X i?", u{T,x) = ^{x), 



(9) 



where u : 



Cu 



Lui 



LUr. 
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with 

Luk(t,x) 
= "^(*'^) + 2 ^i99*)ij{t,x,u{t,x))-^-^{t,x) 

n 
+ X] fi^^^ ^' ^(*' 3;)) Vu(t, x)g{t, X, u{t, x)), 

u(t, j; + /i(t, x, •)) — u(t, a^))-^ — (i, 2;) 

+ / (ufc(t,x + /i(t,x,z)) - Ufc(t,x) - V'/ii(t,x,z)- — {t,x))X{dz), 
Jz ti dxi 

/C = 1, • • • ,771. 

We can assert that 

Theorem 5.1 Assume that the functions f , g, h, F, G and ^ in FBDSDEP 
13) are deterministic and that they satisfy the assumptions (H2)-(H5) and (HI). 
Suppose SPDIE |3j has a unique solution u{t,x) G C^^^{n x [0,T] x M";M"^). 
Then, for any given {t,x), u{t,x) has the following interpretation 

uit,x)=Pt, (10) 

where Pt is determined uniquely by FBDSDEP ^. 
Proof.Applying Proposition 12.41 to u{t,Xt), we obtain 
u{T,XT)-u{t,Xt) 
= j —{s,Xs)ds + J ^f,{s,Xs,Ps,Qs,Ks)-^{s,Xs)ds 

+ j Vu{s,Xs)g{s,Xs,Ps)dWs 

+ {u{s^,Xs_ +h{s^,Xs_,z)) -u{s,Xs_))N(dzds) 

+ 1 -,E(39*hi'^^X.,P.)^^is,X.)ds 

/ {u{s,Xs + h{s,Xs,z)) -u{s,Xs) 
t Jz 

y^hi{^s,Xs,z) — {s,Xs))\(dz)ds. 



Let 



{u{t, Xt),Vu{t, Xt)g{t, Xu u{t, Xt)),u{t, Xt + h{t, Xt, ■)) - u{t, Xt)). 



16 



Because u{t,x) satisfies SPDIE Q, it holds that 

'^{XT)-u{t,Xt) 

= I F{s,Xs,u{s,Xs),Vu{s,Xs)g{s,Xs,u{s,Xs)), 
u{s, Xs + h{s, Xs, •)) - u{s, Xs))ds 

+ / G{s,Xs,u{s,Xs),Vu{s,Xs)g{s,Xs,u{s,Xs)), 
u{s, Xs + h{s, Xs, •)) - u{s, Xs))dBs 

rp 

+ f \/uis,Xs)gis,Xs,uis,Xs))dWs 

+ / (n(s_,X,_ +h{s-,Xs_,z)) -u{s^,Xs_))N{dzds). 

It is easy to check that {u(t,Xt), S/u{t,Xt)g{t,Xt,u{t,Xt)),u{t,Xt + h{t,Xt, •)) — 
u{t,Xt)) coincides with the unique solution to BDSDEP of ([8|). It follows that 

u{t,x) = Pt. 

D 

Remark 5.2 lil(J\) can be called a Feynman-Kac formula for SPDIE ^. Fur- 
thermore, with regard to FBDSDEs driven by Levy process rather than by Poisson 
process, we can get the similar Feynman-Kac formula as lil(J\) for a SPDIE. 

6 Example: a doubly stochastic Hamiltonian system 
with Brownian motions and Poisson process 

Consider the following doubly stochastic Hamiltonian system with Brownian mo- 
tions and Poisson process 

dXt = Hpdt + HgdWt - YtdBt + / HKN{dzdt), 

dPt = -Hxdt-HYdBt + QtdWt+ I Kt{z)N{dzdt), (^^) 

Xo = ^p(Po), PT = '^xiXT), 

where H{X,P,Y,Q,K) : R^ x L2^.)(M) ^ M, $(X) : M ^ R, ^(P) : R ^ M; 
Hp = VpH, <^x = Vx^, ^p = Vp^. The Brownian motions {Wt; 0<t<T} 
and {Bt;0 < t < T} are both assumed to be 1-dimensional. Assume that both 
the derivatives of 2-order of H and the derivatives of 1-order of $ and ^ are 
bounded, H is concave on (P, Q, K) and convex on {X, Y) in the following sense 
/ii > 0, and /X2 > 0: 
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/ —Hxx —Hxp —HxY —HxQ —HxK \ 



Hpx Hpp 


HpY 


HPQ 


HPK 


-Hyx - 


Hyp 


-Hyy 


-Hyq -Hyk 


Hqx Hqp 


Hqy 


Hqq Hqk 


\ Hkx Hkp 


Hky 


Hkq Hkk / 




( l-i-i 





\ 









Ai2 







- 





111 





1 










^2 






^ 





/i2 / 





< 



where V(X, P, Y, Q, K) eR'^x L^^.-, (M), and $ is convex on X: <^xx > 0, VX G R; 
^ is concave on P: ^pp < 0, VP G R. By Theorem 13.21 we claim that this 
doubly stochastic Hamiltonian system (lllh has a unique solution {X, P,Y,Q, K) 

in M2(0,T;R4) x F2,(0,T;]R). 



References 



[lo: 



fii 



fl2 



F. Antonelli, Backward-forward stochastic differential equations, Ann. Appl. Probab. 3 
(1993) 777-793. 

V. Bally, A. Matoussi, Weak solutions for SPDEs and backward doubly stochastic differ- 
ential equations, J. Theoret. Probab. 14 (2001) 125-164. 

G. Barles, R. Buckdahn, E. Pardoux, Backward stochastic differential equations and 
integral-partial differential equations, Stoch. Stoch. Rep. 60 (1997) 57-83. 

J. Cvitanic, J. Ma, Hedging options for a large investor and forward-backward stochastic 
differential equations, Ann. Appl. Probab. 6 (1996) 370-398. 

R. Darling, E. Pardoux, Backward SDE with random terminal time and applications to 
semilmear elliptic PDE, Ann. Probab. 25 (1997) 1135-1159. 

L. Hu, Y. Ren, Stochastic PDIEs with nonlinear Neumann boundary conditions and gener- 
alized backward doubly stochastic differential equations driven by Levy processes, J. Coin- 
put. Appl. Math. 229 (2009) 230-239. 

Y. Hu, S. Peng, Solution of forward-backward stochastic differential equations, Probab. 
Theory Related Fields 103 (1995) 273-283. 

N. Ikeda, S. Watanabe, "Stochastic differential equations and diffusion process", North 
Holland Publishing Company, Amsterdam, 1981. 

J. Ma, P. Protter, J. Yong, Solving forward-backward stochastic differential equations 
explicitly-a four step scheme, Probab. Theory Related Fields 98 (1994) 339-359. 

E. Pardoux, S. Peng, Adapted solution of a backward stochastic differential equation. Sys- 
tems Control Lett. 14 (1990) 55-61. 

E. Pardoux, S. Peng, Backward doubly stochastic differential equations and systems of 
quasilinear parabolic SPDEs, Probab. Theory Related Fields 98 (1994) 209-227. 

E. Pardoux, S. Tang, Forward-backward stochastic differential equations and quasilinear 
parabolic PDEs, Probab. Theory Related Fields 114 (1999) 123-150. 



18 



[13] E. Pardoux, S. Zhang, Generalized BSDEs and nonlinear Neumann boundary value prob- 
lems, Probab. Theory Related Fields 110 (1998) 535-558. 

[14] S. Peng, Probabilistic interpretation for systems of quasilinear parabolic partial differential 
equations, Stoch. Stoch. Rep. 37 (1991) 61-74. 

[15] S. Peng, Z. Wu, Fully coupled forward-backward stochastic differential equations and ap- 
plications to optimal control, SIAM J. Control Optim. 37 (1999) 825-843. 

[16] S. Peng, Y. Shi, Infinite horizon forward-backward stochastic differential equations, Stoch. 
Process. Appl. 85 (2000) 75-92. 

[17] S. Peng, Y. Shi, A type of time-symmetric forward-backward stochastic differential equa- 
tions, C. R. Acad. Sci. Paris, Ser. I 336 (2003) 773-778. 

[18] Y. Ren, A. Lin, L. Hu, Stochastic PDIEs and backward doubly stochastic differential equa- 
tions driven by Levy processes, J. Comput. Appl. Math. 223 (2009) 901-907. 

[19] R. Situ, On solution of backward stochastic differential equations with jumps and applica- 
tions, Stoch. Process. Appl. 66 (1997) 209-236. 

[20] X. Sun, Y. Lu, The property for solutions of the multi-dimensional backward doubly stochas- 
tic differential equations with jumps. Chin. J. Appl. Probab. Stat. 24 (2008) 73-82. 

[21] S. Tang, X. Li, Necessary condition for optional control of stochastic system with random 
jumps, SIAM J. Control Optim. 32 (1994) 1447-1475. 

[22] Z. Wu, Forward-backward stochastic differential equations with Brownian Motion and Pois- 
son Process, Acta Math. Appl. Sinica 15 (1999) 433-443. 

[23] Z. Wu, X. Wang, FBSDE with Poisson process and its application to linear quadratic 
stochastic optimal control problem with random jumps. Acta Automatica Sinica 29 (2003) 
821-826. 

[24] Z. Wu, Z. Yu, Linear quadratic nonzero-sum stochastic diferential games problem with 
random jumps, Appl. Math. Mech. 26 (2005) 1034-1039. 

[25] J. Yin, X. Mao, The adapted solution and comparison theorem for backward stochastic 
differential equations with Poisson jumps and applications, J. Math. Anal. Appl. 346 (2008) 
345-358. 

[26] J. Yin, R. Situ, On solutions of forward-backward stochastic differential equations with 
Poisson jumps, Stoch. Anal. Appl. 21 (2003) 1419-1448. 

[27] J. Yong, Finding adapted solutions of forward-backward stochastic differential equations- 
method of continuation, Probab. Theory Related Fields 107 (1997) 537-572. 

[28] Q. Zhang, H. Zhao, Stationary solutions of SPDEs and infinite horizon BDSDEs, J. Funct. 
Anal. 252 (2007) 171-219. 

[29] Q. Zhu, Y. Shi, X. Gong, Solutions to general forward-backward doubly stochastic differ- 
ential equations, Appl. Math. Mech. 30 (2009) 517-526. 



19 



